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Abstract 

We show that in the quantum query model the complexity of detecting a triangle 
in an undirected graph on n nodes can be done using 0(n 1+ 7 log 2 n) quantum queries. 
The same complexity bound applies for outputting the triangle if there is any. This 
improves upon the earlier bound of 0(n 1+ z). 
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A 1 Introduction 

^ The classical black box complexity of graph properties has made its fame through the noto- 
Rfiously hard evasiveness conjecture of AAnderaa, Karp and Rosenberg (see e.g. [LY]) which 
states that every non-trivial Boolean function on graphs whose value remains invariant under 
^ the permutation of the nodes has deterministic query complexity exactly \ Zj , where n is the 
^ number of nodes of the input graph. 

The question of quantum query complexity of graph properties was first raised in [BCWZ] . 
Here the authors show that the quantum model behaves differently from the deterministic 
model in the zero error case, although an Q(n 2 ) lower bound still holds. 

The bounded error quantum query model is more analogous to the classical randomized 
model. About the latter we know much less, in particular the general lower bounds are 
far from the conjectured f2(n 2 ). For a long time Peter Hajnal's f2(n 4//3 ) bound [H] was the 
best, until it was slightly improved in [CK]. The best general lower bound, fi(n 2 / 3 ), in the 
bounded error quantum case is due to Yao and Santha (not yet published). 

*The research was supported by NSF grant 0105692 
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The quantum query complexity of the Triangle Detection Problem was first treated in 
H. Buhrman, C. Diirr, M. Heiligman, P. Hoyer, F. Magniez, M. Santha and R. de Wolf 
[BDHH], where the authors show that in the case of sparse graphs the folklore n 1,5 upper 
bound can be improved. Their method breaks down when the graph has 8(n 2 ) edges. The 
folklore bound is based on searching the space of all n 3 triangles with Grover's algorithm 
[G]. Since the above search idea is quite natural, one might conjecture that the upper bound 
is the best possible. 

Below we show however, using combinatorial ideas, that a quantum query machine can 
solve the Triangle Detection Problem in time 0(n 1+ ? log 2 n). The lower bound remains Q(n). 

2 Notations 

We denote the set {1, 2, . . . , n} by [n]. A a simple undirected graph is a set of edges G C 

{(a, b) | a,b £ [n]; a ^ &} with the understanding that (a, b) d = (b,a). The complete graph 
on a set v C [n] is denoted by u 2 The neighborhood of a v £ [n] in G is denoted by vq{v). 
Its definition: 

v G {v) = {b\ (v,b)eG}. 
We denote I^gOOI by deg G t>. For sets A, B C [n] let 

G(A,B) = {(a, 6) | a £ A; b £ B; (a, 6) £ G}. 

The following function will play a major role in our proof: We denote the number of 
paths of length two from a £ [n] to b £ [n] in G with t(G, a, b). In formula: 

t(G,a,b) = \ {x | (a,x) £ G; £ 

is the number of triangles in G. We define an operation on graphs. For a graph G C [n] 2 
define 

G<*> = {(a, 6) £ [n] 2 | t(G,a,b) < t}. 

3 The Triangle Detection Problem (TDP) 

Input: An undirected graph G with vertex set [n], given by its incidence matrix, G(a,b). 
(G(a, b) is a symmetric zero-one valued matrix with all zeros in the diagonal.) 
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Output: "No," if G is triangle free. Else, a set {a, b, c} such that (a, b), (b, c), (c, a) G C 

Any classical machine that solves the Triangle Detection Problem (TDP) has to query 
all edges of the graph, and even a classical randomized machine has to query a constant 
fraction of all edges. In contrast, it has been known that a quantum query machine can 
successfully output a triangle of G, if exists, making only 0(n L5 ) quantum queries to G, 
with high probability. In this section we give an algorithm that improves on the exponent 
1.5. 

The basis states of a quantum query machine for the TDP are of the form \a, b) |c), where 
1 < a < b < n and c G S is specific to the the machine. The starting state is |1, l)|c ) for 
some distinguished c G S. A query step is a unitary operator O = Og acting on the basis 
elements as: 

G \a,b)\c) = (-lf^\a,b)\c). 

The computation is a sequence U Og ■ ■ ■ U t -iOaU t of query steps and non-query steps. 
Each non-query step is an arbitrary unitary operator on the state space with the only restric- 
tion that it does not depend on G. The output is read from a measurement of the final state, 
and needs to be correct with probability at least 1 — e. We will describe our algorithm as 
a one that runs on a classical machine that calls quantum subroutines, but makes classical 
queries as well. Although this type of the machine seems more general than the one we 
explained above, the two models can be shown to have the same power. 

The quantum query complexity of the Triangle Detection Problem as well as of many of 
its kins with small one-sided certificate size are notoriously hard to analyze, because one of 
the main lower bounding methods breaks down near the square root of the instance size: 

Lemma 1 // the 1-certificate size of a Boolean function on n Boolean variables is k then 
the general weighted version of the quantum adversary method ( also known as the method of 
Ambainis) in the sense of [BSS] can prove only a lower bound of (l — 2\Je(l — e)j \fnk. 

Problems with small certificate complexity include various collision type problems such 
as the 2-1 collision problem and the element distinctness problem. Using the degree bound of 
R. Beals, H. Buhrman, R. Cleve, M. Mosca and Ronald de Wolf [BBCM] for the 2-1 collision 
problem the first polynomial lower bound was shown by Aaronson [Aa]. Shi [S] showed tight 
H(n 1//3 ) lower bound and Kutin [K] removed the "small range assumption". For element 
distinctness, an Vt{n l l 2 ) lower bound follows by an easy reduction from Grover's search. 
Shi got fi(n 2 / 3 ) and Ambainis [A] removed the "small range assumption." With a recent 
ingenious algorithm of Ambainis [A2] the complexity of the element distinctness problem 
was finally determined to be in 6{n 2 ^). In a sequel of this paper using Ambainis's new 
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technique with Santha and Magniez [MSS] we could improve the present upper bound for 
the TDP to n 13 . The algorithm presented here is based on three combinatorial observations: 
The first observation is due to H. Buhrman, C. Diirr, M. Heiligman, P. Hoyer, F. Magniez, 
M. Santha and R. de Wolf [BDHH]: 

Lemma 2 If graph G\ is known then detecting a triangle with at least one edge G\ can be 



Proof: The proof of the above statement is exactly the same as that of Theorem 8 in 
[BDHH], and is based on the Amplitude Amplification technique of G. Brassard, P. Hoyer, 
M. Mosca, A. Tapp [BHMT]. □ 

Note that the lemma can also be proven using Grover search, and we loose only a factor 
of log n in the analysis. Throughout the paper we do not try to optimize log n factors, since 
we conjecture that the exponent we present here is not tight. Perhaps the most crucial 
observation to the algorithm is the following trivial-looking one: 

Lemma 3 For every v G [n], using O(nlogn) queries, we either find a triangle in G or 
verify that G C [n] 2 \ vq{v ) 2 with probability 1 — -K,. 

Proof: We query all edges incident to v classically using n — 1 queries. This determines 
vg(v)- With Grover's search [G] we find an edge of G in v G {v ) 2 , if there is any. In order to 
achieve 1 — ^3 success probability we use the following safe variant of the search: 

Definition 1 (Safe Grover Search) Let c > 1 be an arbitrary constant. The Safe Grover 
Search on a database of N items is the usual Grover search iterated O(logiV) times indepen- 
dently. The safe Grover search uses O c (N log N) queries and it finds the (type of) item we 
are looking for, if it is the database, with probability at least 1 — \. 



This lemma with the observation that hard instances have to be dense, already enable 
us to show that the quantum query complexity of the TDP is o(n L5 ), using the lemma of 
Szemeredi. However another fairly simple observation can help us to decrease the exponent: 

Lemma 4 Let k = |~4n e logn] 7 and let v±, V2, ■ ■ ■ , Vk randomly chosen from [n] (with no 
repetitions). Let G' = [n] 2 \ U k i=l v G {v) 2 . Then 



determined with 0(Jn\G fl G'\) queries. 



□ 
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Proof: Let us first remind the reader about the following lemma that is useful in many 
applications: 



Lemma 5 Let X be a fixed subset of [ra] of size pn and Y be a random subset of [raj of size 
qn, where p + q < 1. Then the probability that X fl Y is empty is (1 — pq) n ( l± °(P +<?' . 

Proof: The probability we are looking for is estimated using the Stirling formula as 
CV) _ [ra(l-p)]!H![ra(l- g )]! 

n 
nq 



;i-p)(i-g) 
\ (i-P-?) 



[rag]![ra(l — p — q)]\n\ 

1 _ p )l-p(l _ g)l-<? 



(1 -p - g) 1 "^ 



l±o(l)) = (l- p? )"( 1±0 (f 3 +« 3 +Vn))_ 



□ 

Consider now a fixed edge (a, 6) such that t(G, a, b) > ra 1_e . The probability that (a, 6) G 
G" is the same as the probability that the set X = {x e [ra] | (x, a) e G A (x, 6) G G} is 
disjoint from the random set {v±, t> 2 , . . . , ffc}. Notice that \X\ = t(G, a, b). By Lemma 5 we 
can estimate now this probability as 

/ _ toMog^ x ^-.\"C«C» = / _ 4 log n \ ""-""» < e _ 3l0g „ = n _ 3 ^ 
1 ra ra / 1 ra / 



Then the lemma follows from the union bound, since the number of possible edges (a, b) is 
at most ra 2 .D 

Our algorithm will have three parameters e = |, e' = 5 = j, and it either returns a 
triangle during its run or it: 

1. Randomly creates a graph G' such that GCG'C G^ 1 ^ with high probability; 

2. In a number of steps classifies the edges of G' into T and E such that T contains only 
a small number of triangles and E fl G is small; 

3. Searches for a triangle in G among all triangles inside T; 

4. Searches for a triangle of G intersecting with E; 

5. In case steps 3 and 4 run without success the algorithm returns "No." 
See the algorithm in details in Figure 1. 
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Command 



Complexity 



1. Set k = |~4n e logra] Select {vi, v 2 , ■ ■ ■ , v^} randomly from [n] and 0(n l+e logn) 
query all edges of the form (vi, v) for 1 < % < k and v G [n]. 

2. For each Vi check with the Safe Grover search if G n voipi) 2 is 0(n 1+e (logn) 2 ) 
empty. If some 1 < % < k it is not, after verification return 

(v,w), (vi,v), (vi,w), where all three edges are in G. 
Otherwise define G' = [n] 2 \ U k i=1 u G {vi) 2 . 

3. Set T = 0, E = 0. 

4. Find an edge (v, w) G G' for which t(G', v, w) < [n 1_e '] if exists, 
and put this edge to T, delete it from G' . Keep repeating this 

until no such edge is found. 

5. Pick a vertex v of G' with non-zero degree and choose one of 0(n s logn) 
the 1. low degree hypothesis: \vq{v)\ < 10[n 1_<5 ]; 2. high degree 
hypothesis: \ug(v)\ > 0.1 [n 1_<5 ] ; such that the chosen hypothesis 

holds with probability at least 1 — ^3- 

6. If in 5. we chose the low degree hypothesis, add all edges of G' 
with end point v into E and delete them from G'. 

7. If in 5. we chose the high degree hypothesis, then find all vertices 0(n log n) 

def 

in A = vg{v). First run the Safe Grover to determine if there 
is an edge in vg{ v ) 2 - If one is found, after verification, output 
the triangle of G induced by v and this edge. If no triangle is 
found add all edges in G'(A, A') to E, where A' = u G >(v). Also, 
delete these edges from G'. 

8. Repeat 4.-7. until G' becomes empty. See analysis 

9. Try to find a triangle inside T using Grover search, verify and V n 3 ~ 6 ' log n 
output if found. If not, go to 10. 

10. Search for a triangle in G with non-empty intersection with E, \Jn\G fl E\ 
verify and output if found. 

Figure 1: Quantum query algorithm for the Triangle Detection Problem 
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4 Analysis 



Lemma 6 The algorithm for the TDP described in the previous section runs using 

0(n 1+e (\ogn) 2 + n 1+s+t ' logn + Vn 3 " 6 ' log n + v / n 3 - min ( <5 ' £ - <5 - e ')). 

quantum queries and returns "No" with probability one if there is no triangle in G, otherwise 
returns a triangle of G with some constant positive probability. 

Proof: 

We prove the correctness and the bound on the running time together. Clearly, if there is 
no triangle in the graph, the algorithm outputs "No," since the algorithm places a triplet into 
the output register only after checking that it is a triangle in G. Therefore the correctness 
proof requires only to calculate the probability with which the algorithm outputs a triangle 
if there is any. 

Step 1: 

For each t>j we use n — 1 deterministic queries to find VQ{vj). Since 1 G [k], where 
k = 0(n e \ogn), Step 1 costs 0(n 1+e log n) queries. 

Step 2: 

For each Vi we use Safe Grover to find out if ^g(^) 2 is empty. A single run of Safe Grover 
takes O(nlogn) queries, thus the total number of queries made is 0(n 1+e (log n) 2 ). We now 
claim that with probability \ — 0{\) we have that G C G' C G^ n ~*\ or else a triangle of G 
is found. 

Indeed, if no triangle is found then the Safe Grover for each Vi with probability 1 — \ 
verifies that Gnz/ G (^) 2 = 0. So G' D G holds. That G' C G^~^ holds follows from Lemma 
4. 

Step 4-' 

Realize that after steps 1. and 2. we know G f , so Step 4 costs us no queries. One of 
the important, but simple observations is that even repeated applications of step 4 creates 
a graph T with a small number of triangles. This is shown in the following lemma: 

Lemma 7 Let H be a graph on [n] . Assume that a graph T is built by a process that starts 
with an empty set, and at every step either discards some edges from H or adds an edge 
(a, b) of H to T for which t(H, a,b) < t holds. For the T created by the end of the process 
we have t(T) < (™)r. 
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Proof: Let us denote by T[i\ the edge of T that T acquired when it was incremented for the 
2 th time, and let us use the notation H l for the current version of H before the very moment 
when T[i] = (a h h) was copied into T. Since {T[i], T[i + 1] . . .} d = T l C H\ we have 

t(T\ ai ,h) < t(H\a h h) < r. 

Now the lemma follows from 



t(T) = Y,t(T\aM<(^jT, 



since i can go up to at most (™) • 1=1 
Step 5: 

In this step we use an obvious sampling strategy: Set a counter C to 0. Query \n s ~\ 
random edge candidates from v x [n] (the constants are important). If there is an edge of G 
among them, add one to C. Repeat this process K = c logn times, where c is a sufficiently 
large constant. Accept the low degree hypothesis if by the end C < K/2, otherwise accept 
the large degree hypothesis. Then 

1. The probability that deg G (i>) > lOn 1 ^ 5 and the low degree hypothesis is accepted is at 
most 1/n 3 . 

2. The probability that deg G (v) < O.ln 1 ^ 5 and the high degree hypothesis is accepted is 
at most 1/n 3 . 

Indeed, using Lemma 5, considering a single round of sampling the probability that our 
sample set does not contain an edge from G even though deg G (v) > lOn 1 ^ 5 is 

/ lOn 1 - 5 n *x"(i+*(i)) / ion 

1 x— = 1 <0.1 

\ n n J \ n / 

Similarly, the probability that our sample set contains an edge from G even though deg G (i>) < 
O.ln 1 - 6 is 

/ O.ln 1 ' 5 n 5 \ n(1+o{1)) ( 1 V l(1+o(1)) 

1-1 x— =1-1 < 0.2. 

\ n n ) \ 10n/ 

Now for K = Co logn rounds, where Co is large enough the Chernoff bound gives the 
claim. 
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Step 6: 

Since we know G' the operation costs us no queries. 
Step 7: 

Finding out u G (y) costs us n — 1 classical queries. Finding out if z/ G (t>) 2 n G = costs us 
O(nlogn) queries, using the safe Grover. 

Step 8: 

The key to estimating the complexity of this step is an upper estimate on the number of 
executions of Step 7. In turn, this is done by lower bounding \G'(A, A')\. For each xeAwe 
have t(G',v,x) > rf \ otherwise in Step 4. we would have classified (v,x) into T. A triangle 
(v, x, y) contributing to t(G', v, x) contributes with the edge (x, y) to G'(A, A'). Two triangles 
(v, x, y)and (v, x', y') can give the same edge in G'(A, A') only if x — y' and y = x'. Thus: 

\G'(A,A')\>± t(G?,v,x)>\A\T?-*/2. (1) 

x£uq{v) 

Since we executed Step 7 only under the large degree hypothesis on v, if the hypothesis 
is correct, the right hand side of Equation 1 is at least 0.1n 1 ~ s n 1 ~ t /2 = £l(n 2 ~ s ~ € ). Since 
G' has at most (™) edges, it can execute Step 7 at most 0(n 5+e ') times. 

How about the number of executions of Step 5? Every execution of Step 5 leads to either 
the execution of Step 6 or that of Step 7. We have already seen that the number of executions 
of Step 7 is in n s+t , which is easily in 0{n) (and it would not help to choose the parameters 
otherwise). We claim that each vertex is processed in Step 6 at most once. Indeed, if a 
vertex v gets into Step 6, its incident edges are all removed, and its degree in G' becomes 
making it ineligible for being processed in Step 5 again. Thus the total complexity of Step 
8 is: 

n 5+t O (n log n) + 0(n)0(n 5 log n) . 

Step 9: 

By Lemma 7 we have that there are at most n 2 n 1_e ' triangles are in T. T is a graph 
that is known to us, and so we can find out if one of these triangles belong to G in time 
0(Vn 3 ~ e ' logn), using Safe Grover. 

Step 10: 

In order to estimate G D E observe that we added edges to E only in Steps 6. and 7. 
Assume the Safe Grover worked throughout the whole algorithm. In each execution of Step 
6. we added at most lOn 1 ^ 5 edges to E, and we had 0(n) such executions that give a total 
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of 0(n 2 ~ 5 ) edges. The number of executions of Step 7 is 0(n 5+e '). Our task is now to bound 
now the the number of edges of G each such execution adds to E. 

We estimate \G n G'(A,A')\ from the A' side. This is the only place where we use the 
fact that G' C G^ n e K For every x G A' we have t(G, v, x) < n 1 ^ t . On the other hand every 
edge (y, x) y G A, x G A', (y, x) G G' creates a (i>, x)-based triangle. Thus 

\Gr)G'(A,A')\ < lA'ln 1 - 6 < n 2 ~ e . 

Therefore the total number of edges of G Step 7. contributes to E is n 2 ~ e+s+e ' . In 
conclusion, 

\GDE\ < 0(n 2 ~ 5 + n 2 - e+5+e '). 
By Lemma 2 the complexity of finding a triangle in G that contains an edge from E is 

O (y n 3-min(S,e-S-e>)^ _ 

4.1 Conclusions of the Analysis 

The probability that not all the Safe Grovers run correctly can be upper bounded by ^, for 
large n, using the union bound. Thus in the sequel we assume that none of the Safe Grovers 
fails. From the analysis we conclude that the total number of queries is upper bounded by: 

0{n l+e logn + n 1+e (logn) 2 + (n 1+5+t ' logn + n 1+s logra) + Vn 3 ~ e ' \ogn + y/ n ^m(8,e~5^')y 

With e = |, e' — S — | this gives 0(n 1+ ^ log 2 n) for the total number of queries. When 
the Safe Grover does not fail, the structural assumptions we made based on it hold. G' 
eventually lends all its edges to T and E. Since G C G', every triangle in G either has to 
be contained totally in T or it has to have a non-empty intersection with E. Finally, we 
mention that in case one of the Safe Grovers goes wrong, the algorithm cannot run beyond 
the given time bound: quantum query machines by definition execute a prescribed number 
of steps. 

Acknowledgments: The author thanks Frederic Magniez and Miklos Santha for bringing the 
problem and Lemma 2 to his attention and Andris Ambainis for references. 

References 

[Aa] Scott Aaronson: Quantum lower bound for the collision problem. STOC 2002: 
635-642 



10 



[A] Ambainis: Quantum Lower Bounds for Collision and Element Distinctness with 
Small Range, quant-ph/0305179. 

[A2] Ambainis: lecture on the complexity of collision finding (communicated by Frederic 
Magniez) 

[BSS] Howard Barnum and Michael Saks and Mario Szegedy: Quantum decision trees 
and semidefinite programming. CCCC 2003 

[BBCM] Robert Beals, Harry Buhrman, Richard Cleve, Michele Mosca, Ronald de Wolf: 
Quantum Lower Bounds by Polynomials. FOCS 1998: 352-361 

[BCWZ] Harry Buhrman, Richard Cleve, Ronald de Wolf, Christof Zalka: Bounds for Small- 
Error and Zero-Error Quantum Algorithms. FOCS 1999: 358-368 

[BHMT] Gilles Brassard, Peter Hoyer, Michele Mosca, Alain Tapp: Quantum Amplitude 
Amplification and Estimation, quant-ph/0005055 

[BDHH] Harry Buhrman, Christoph Diirr, Mark Heiligman, Peter Hoyer, Frederic Magniez, 
Miklos Santha, Ronald de Wolf: Quantum Algorithms for Element Distinctness. 
IEEE Conference on Computational Complexity 2001: 131-137 

[CK] Amit Chakrabarti, Subhash Khot: Improved Lower Bounds on the Randomized 
Complexity of Graph Properties. ICALP 2001: 285-296 

[G] Lov K. Grover: A Fast Quantum Mechanical Algorithm for Database Search. STOC 
1996: 212-219 

[H] Peter Hajnal: An n 4//3 lower bound on the randomized complexity of graph prop- 
erties, Combinatorica, 11(1991), 131-143 

[LY] L. Lovasz and N. Young. Lecture notes on evasiveness of graph properties. Technical 
Report CS-TR-317-91, Computer Science Department, Princeton University, 1991. 

[K] Samuel Kutin A quantum lower bound for the collision problem quant-ph/0304162 

[MSS] F. Magniez, M. Santha, and M. Szegedy: An 0(nl.3) quantum algorithm for the 
triangle problem. Technical Report quant-ph/0310134, arXiv, 2003. 

[S] Yaoyun Shi: Quantum Lower Bounds for the Collision and the Element Distinctness 

Problems. FOCS 2002: 513-519 



11 



5 Appendix 



Proof of Lemma 1: 



For a matrix M, let A(M) denote the greatest eigenvalue of M. Recall the general form 
of Ambainis from [BSS]: 

Theorem 1 Let f : S — > T be a partial boolean function with S C {0, 1}". Let T be an 
arbitrary SxS nonnegative symmetric matrix that satisfies F[x, y] — whenever f (x) = f(y). 
For i G {!,...,«} let Tj be the matrix: 



Fi[x,y] = 

Then: 



0ifxi = yi 
Y[x,y\ ifxi ^ yi. 



(l - 2Je(l-e)) A(r) 

Above QQC e (/) denotes the e error quantum query complexity of /. We need to work 
in the case when SCE" and T = {0, 1}. Let A = Let k be the 1-certificate size of 

/. For an input x G A let A x be its 1-certificate. Let 



V[x,y\' 



T[x,y] if f(x) = 1 and f(y) = 
otherwise. 



Let v E H D be the unit vector with A(r) = vTv*. Notice that v is positive. Because 
r[x, y] is whenever f(x) = f(y) and T is symmetrical, it holds that vT'v* = A/2. 

For 1 < i < n let V{ be the vector which equals to v on those coordinates x where i G A x 
and otherwise. We now have: 



(vi,v) = (vi,Vi) for 1 < i < n. (2) 

vi + . . . + v n < kv; (3) 

v^v + v 2 T 2 v + ... + v n F n v > vT'v = A/2; (4) 

A(r) 



max!<j< n A(rj) 



< 2(K| + ... + K|). (5) 
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Here the first inequality is meant entry-wise, and Inequality (5) follows from Inequality 
(4) and from ViTiV = \vi\-^TiV < \vi\\(Ti). Now from (2) and (3) it follows that 



EN < 



nJ2\ v i\ 2 



, n^2(vi\v) < yn(kv\v) = vnk. 
\ i=i 



The above together with (5) gives the lemma. □ 
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